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Abstract 

The Cosmic Microwave Background can provide information regard¬ 
ing physics of the very early universe, more specifically, of the matter- 
radiation distribution of the inflationary era. Starting from the effective 
field theory of inflation, we use the Goldstone action to calculate the three 
point correlation function for the Goldstone field, whose results can be di¬ 
rectly applied to the field describing the curvature perturbations around 
a de Sitter solution for the inflationary era. We then use the data from 
the recent Planck mission for the parameters and which 

parametrize the size and shape of non-Gaussianities generated in single 
field models of inflation. Using these known values, we calculate the pa¬ 
rameters relevant to our analysis, /vn, ^^e speed of sound Ca 

which parametrize the non-Gaussianities arising from two different kinds 
of generalized interactions of the scalar field in question. 


1 Introduction 

The Cosmic Microwave Background(CMB) is an important probe of physics 
of the early universe. We study the phenomenon of primordial inflation and 
its observational signatures on the CMB spectrum(for review see [IllHIS]). In 
standard inflationary scenarios, the universe should be very close to a Gaussian 
random field. The two point correlation function and its Fourier transform, the 
angular power spectrum for most models of inflation give similar prediction of a 
scale invariant, adiabatic, Gaussian spectrum. Therefore, in order to distinguish 
between the competing models of inflation and to constrain the parameters 
common to these models, we look at the bispectrum predictions from these 
models, and deviations from Gaussianity. Such a calculation of three-point 
function was first done by Maldacena in [4] . The theory of single field slow roll 
inflation achieves accelerated expansion by means of a scalar field slowly rolling 
down a potential. It predicts non-Gaussianities that should be too small for 
observation. However, other models of inflation predict larger departures from 
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the Gaussian spectrum, and have their characteristic types of non-Gaussianities. 
Any observation of large non-Gaussianity from the Planck mission thus allows 
us to constrain these models. Moreover, the angular bispectrum has different 
shapes, and the different models of inflation show peaks for distinct shapes. 
As elaborated in B IS] and references therein, various models of primordial 
non- Gaussianity are known as local, equilateral, orthogonal or folded models in 
literature. Different aspects of physics of the early universe appear in different 
shapes of the three point function. 

• Local Non-Gaussianity appears in multi-field models of 
inflation due to interactions which operate on superhorizon scales. 

• Equilateral Non-Gaussianity includes single field models with non- 
canonical kinetic term such as k-inflation or Dirac-Born- Infeld inflation 
models [n m nn [n] characterized by more general higher derivative 
interactions of the inflaton field such as ghost inflation and models arising 
from effective field theories. 

• Folded Non-Gaussianity include single field models with non Bunch 
Davies vacuum. 

• Orthogonal Non-Gaussianity may be generated in single field mod¬ 
els of inflation with a non-canonical kinetic term or with general higher 
derivative interactions. The orthogonal form is constructed in such a way 
that it is nearly orthogonal to both local and equilateral forms. [7] 

The plan of this article is as follows: In the next section we review the 
effective field theory model of inflation. After that we explicitly calculate 
2-point and 3-point function in section 3 and finally conclude in section 4. 

2 The Effective Field Theory of inflation 

Gheung, Greminilli, Fitzpatrick, Kaplan and Senatore have used the effective 
field theory approach in [S] to describe the theory of fluctuations around an 
inflating cosmological background. While the inflaton field (/> is a scalar under 
all diffeomorphisms, the perturbation S(j) is a scalar only under spatial diffeo- 
morphisms, and transforms non-linearly with respect to time diffeomorphisms, 

t ^ t + x) 

Sc/> S(f> + 

We can describe the perturbations during inflation directly around the time- 
evolving vacuum where the time diffeomorphisms are non-linearly realised. In 
unitary gauge, the most generic Lagrangian with broken time diffeomorphisms 
and unbroken spatial diffeomorphisms around a flat FRW with Hubble param¬ 
eter ff(t) is given by [8], 

S = J + JI) + + 1)" 
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{g^^ + iy + ...+ 


Miit) 


iSK^r 


As elaborated by Senatore, Smith and Zaldariagga, in inflation there is a 
physical clock that controls the end of inflation, so that time translations are 
spontaneously broken, and there is a Goldstone boson associated with the sym¬ 
metry breaking. The Lagrangian of the Goldstone boson is highly constrained 
by the symmetries of the problem, in this case the fact that spacetime is ap¬ 
proximately de Sitter, and -^ << 1. The Goldstone boson, tt can be thought 
of as being equivalent, in standard models of inflation driven by a scalar field, 
to the perturbations in the scalar field Scj)- The relation valid at linear order 
is TT = S(j)/(j), where (j) is the speed of the background solution. The Gold¬ 
stone boson is related to the standard curvature perturbation ^ by the relation, 
^ = —Htt, which is valid at linear order and leading order in the generalised 
slow roll parameters. The most general Lagrangian for the Goldstone boson is 
given by imm] The general action for the perturbation tt is 





(9i7r) 


M^iH ( 1-1 


— fr 


{diTif 


( 1 ) 


To arrive at the Lagrangian for the Goldstone mode from the most generic 
Lagrangian in unitary gauge, one performs a time diffeomorphism with param¬ 
eter and promotes the parameter to a field, '!r(t,x), which shifts under 

time diffeomorphisms, Tr{t,x) —>■ Tr{t.x) + ^°{t,x). This scalar, tt, is the Gold¬ 
stone mode which non-linearly realises time diffeomorphisms and describes the 
scalar perturbations around the FRW solution [8]. 


3 Calculation of Non-Gaussianities 

We use the action described above in Eq. [T] here to calculate the bispectrum. 

Solution of quadratic action Considering only upto the quadratic terms 
the action is. 


52 = / d^xi-M^i)H 




- a{diTr)‘ 


The Euler- Lagrange equation for the field tt becomes, 
d 


dt 


TT 1 - 2M^iHaV^n = 0 


Decomposing tt into momentum modes using Fourier transform, 

7r(r, x) = [ d^kTTk{T)e^^’^ 
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we get 


d I . 

dt c2 


7r'k I + 2MpiHak TT\^ = 0 


Now defining the Mukhanov-Sasaki variable 

Vk = ZT^k 

where 


z = 


a\/-2H M, 


pi 


and transforming everything to conformal time r, such that dt = a dr we get 


< + ( k^d - 


Wk = 0 


where ' represents derivative with respect to conformal time. Now the strongest 
time dependence is contributed by scale factor a so if we take H and Cs to be 
varying slowly, we can write ^ = dr- perfect de Sitter space a = — 


so 


So the equation reduces to 


k d “ = 0 


( 2 ) 


which has a solution 


For perturbation tt, 


^ 

\/2kcs Cgkr 


= id + ikcsT) ^-ikc,T 
z 2y/e-Csk k Mpi 


( 3 ) 


where e = —is the usual slow roll parameter and we have used a = — 


Differentiating it with respect to r we get 

/ i 


2y/ekcs Mpi 


{ki 


J2^^—ikcsT\ 


( 4 ) 


Now we need to think of cubic terms in the action as the perturbation to 
quadratic action. In other words cubic terms will contribute to interaction 
part of the Hamiltonian. Upto the cubic order we have 'Hint = ~ k^int so 
the interaction Hamiltonian in leading order becomes 


H^ntit) = - d^a; d 


Milk 


1 - 


•d — TT 


. {diTrf 


- -Mid 


( 5 ) 


This has two terms: one is dd and the other one is a'k{diTTd ■ We will calculate 
the correlation function corresponding to each separately. 


4 














Calculation for a^TT^ term We can expand tt into creation and annihilation 
operators as 

TT = TTkOk + TT^aj^ 

Weinberg’s in-in formalism m gives us a very nice way to calculate correlation 
functions. According to it any correlation function ((bb(t))) is given by 

pt ptN /‘t2 


If we consider only upto first order in interaction we get 
{W{t))=if dt{[Hintit),W{t)]) 

Jto 


We will take to to be the beginning of the universe which will become — oo 
once we convert everything to conformal time. We do the integral upto r = 0. 

(7r(ki,0)7r(k2,0)7r(k3,0)) =-z / adr 

J —OO 


d'’’a;([7r(ki, 0)7r(k2, 0)7r(k3,0), 7r"(T, x)]) 


( 6 ) 


Now again decomposing tt' into Fourier modes 


(7r(ki,0)7r(k2,0)7r(k3,0)) =-z(27r)^y adr J dk 4 J dk^ J dke{[Tr(ki,0)7r{k2,0)Tr{'k3,0), 

7r'(k4,T)7r'(k5,r)7r'(k6,r)])5^(k4 + ks -h ke) 

Now using Wick contraction and considering only connected diagrams 

M 

(7r(ki,0)7r(k2,0)7r(k3,0)) = -6z(27r)^5^(ki+k2-fk3) / adr(7rki (0))(7rk^( t))* 

J —OO 

7''k2(0))(7rkJr))*(7rk3(0))«3(T))* -f C.C. 

Putting values of TTk and 7r(^ from Eq. |3] and SI we get 

3zc^ 


(7r(ki,0)7r(k2,0)7r(k3,0)) = - 


32e3fcifc2A:3M6 


(27r)353(ki-p ka -b kg) 


pi 


/ ar3 

^i{ki-\-k2-\-k3)csr _ ^—i{ki+k2-\-k3)csT 

y —OO 

■ 


dr 


or putting a = —-4- we get 


(7r(ki,0)7r(k2,0)7r(k3,0)) = -b 


3^3 


32e^Hkik2k3M^i 




^0 ^0 

^ 2 ^i{ki+k 2 +k 3 )ceT^^ _ / ^ 2 ^i(ki+k 2 + k 3 )csT^^ 
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Now this integral doesn’t converge. But we can get rid of divergences using 
a slightly imaginary time axis while doing integration as follows: 


(7r(ki,0)7r(k2,0)7r(k3,0)) = + 


32e^Hkik2k3M^i 

^2gi(fci+fc2 + fc3)csr^^ _ / 

J — oo(l+ie) 




(2^)^d^(k; 


kg) 


^-i{kl+k2 + k3)CeT 


' —oo(l—ie) 

which gives 

(^(ki,0K(k2,0)7r(k3,0)) = - St^Hk^klk.K^M^ (2^)^d"(ki+k2+k3) (7) 

pi 

where K = ki k 2 k^. 


Calculation for atr(9i7r)^ term Now doing this calculation for a‘k{diTT)'^ term, 
we get 

(7r(ki,0)7r(k2,0)7r(k3,0)) =-i f adr / fi^x([7r(ki, 0)7r(k2,0)7r(k3, 0), 7r'(ai7r)^(r, a;)]) (8) 

J —CO J 

(7r(ki,0)7r(k2,0)7r(k3,0)) =-i(27r)^y adr J ( 1^4 J y dk6([7r(ki, 0)7r(k2,0)7r(k3, 0), 

7r'(k4, T)9i7r(k5, r)5i7r(k6, r)]) d^(k4 + ks + ke) 

Now using Wick contraction and considering only connected diagrams 

(7r(ki,0)7r(k2,0)7r(k3,0)) = i(27r)^(5^(ki+k 2 d-ka) / adr 

J —OO 

[2k2.k3(7rki (0))(7rk^ (0))(7''k2 (t))* (7rk3 (0))(7rk3 (r))* 

+2ki.k3(7rki (0))(7rki (r))*7rk3 (0))«^ (t))*( 7rk3 (0))(7rk3 (r))* 

+2ki .k2(7rki (0))(7rk3 (T))*7rk3 (0))(7rk3 (r))* (rka (0))(7rl,^ (r))*] 

+C.C. 

Now since ki, k 2 and ka form the sides of a triangle, 2k2.k3 = {kf — — k^) 

and similarly for other two permutations. Now we put the values of TTk and 
from Eq. [3]and|4l and integrate. After some algebra we get, 

(7r(ki, 0)7r(k2, 0)7r(k3, 0)) = " oo a (27r)^(5^(ki + k 2 + ka) 


where 


32e3c2i7(M)3M6 
+ 12M^ -4KLM-4K^L^ + IIM - 3LK^ 

^3 

K = ki + k2 + k3 
L = kik2 + ^ 2^3 + k^ki 
M = kik2k3 


( 9 ) 
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( 10 ) 

( 11 ) 

( 12 ) 








Total Bispectrum: We can get total bispectrum by adding the above two 
terms multiplied by the coefficients in the original action.This gives: 


(7r(A:i,0)7r(fc2,0)7r(A:3,0)) = -(2^f <5"(ki + ka + k3)[- ( ( 1 - ^ 


1 


+ M^iH ( 1 - - 


ciJ 32e3c2iJ(M)3M6 


c^J 3 V 8e3i7fcifc2fc3^3M® 

+ 12M2 -4:KLM-4:K‘^L‘^ + IIM - 3LK'^ 
^3 


In the effective field theory model this tt can be related to r that remains 
unchanged after horizon crossing as 


So, 


or 


C(t,x) = -iJ7r(f,x) 

(C(A:i,0)C(A:2,0)C(fc3,0)) = -H^{Tr{ki,0)Tr{k2,0)'!r{k3,0)) 


(C(A:i,0)C(fc2,0)C(A:3,0)) = (27r)353(-]^^^j^^_l_]^g^ [^* 3 (^ 1 , fe, ^ 3 ) +-FV(ai77)2(fci, ^2, ^ 3 )] 
where 


F^a(fci, fc 2 , fc3) = - ( MliH 


1 


3iJ2 


c2; 3 y 8e^kik2k3K^M% 


(13) 


pZ 


and 


^iT( 0 i 7 r )2 (^ 1 ; ^2 7 ^ 3 ) (1 2 ) ^ 


H^H 


cV 32e3c2(M)3M4 
fK^ + 12M2 -4:KLM-4:K‘^L‘^ + IIM - iLK^ 


V 




(14) 


Power Spectrum : Let’s now define power spectrum as 

Ac(fc) 


(C(ki,0)C(k2,0)) = (27r)3<l3(ki-k2) 


fc3 


Also 


or 


So 


(C(ki, 0)C(k2, 0)) = i/2(2^)3^3(j^^ _ |2| 

1 


(C(ki,0)C(k2,0))=iL^(27r)353(ki-k2) 


= 4efc3c,M2 


Aek^CsM^i 
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Experimental constraints : We define /nl as 


Inl 


5 F{k, k, k) 
18 A^(/c)2 


For Tr(9i7r)^ term 


and for term 


JNL 


85 



1 






Now this shape doesn’t match equilateral or orthogonal shape. It is a combina¬ 
tion of both. The exact relation between and have 

been given in [12] as 


nequil. 

JNL 

porthoq. 

Inl 


1.040 1.210 \ 

-0.03951 -0.1757 ) I 


Inverting these relations, we get 

= 1 . 3022 /^«“'' + 8 . 9682 /°’-^°®- 
/*! = -0.2928/X'' - 7.7082/°"^°®' 

Planck [TT] sets the constraints as —86 < < 54 and —67 < < —I 

( 68 % CL statistical) which implies constraints as 

- 712 < < 61 (15) 


- 8.103 < < 541 

This gives for speed of sound 


> 0.019 


(16) 

(17) 


4 Conclusion 

We use the effective field theory model of inflation and analytically calculate 
the 2- point and the 3-point functions for the field perturbations. Then we use 
the latest Planck data m to constrain various parameters of the model. 

This similar calculation can be found in |9] where authors have used WMAP 
data to constrain the parameters. Our final form of three point function matches 
with that. Moreover here we show all the calculation steps explicitly and use 
the latest Planck data. Also the constraint on speed of sound Cg was found by 
Planck collaboration in m as Cs > 0.021 for DBI model. 
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